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ABSTRACT 



Context. Rings or annulus-like features have been observed in most imaged debris discs. Outside the main ring, while some systems 



(e.g.,/3 Pictoris and AU Mic) exhibit smooth surface brightness profiles (SB) that fall off roughly as 



, others (e.g. HR 4796A 



and HD 139664) display large drops in luminosity at the ring's outer edge and steeper radial luminosity profiles. 

Aims. We seek to understand this diversity of outer edge profiles under the "natural" collisional evolution of the system, without 

invoking external agents such as planets or gas. 

Methods. We use a multi-annulus statistical code to follow the evolution of a collisional population, ranging in size from dust grains 
to planetesimals and initially confined within a belt (the "birth ring"). The crucial efi'ect of radiation pressure on the dynamics and 
spatial distribution of the smallest grains is taken into account. We explore the dependence of the resulting disc surface brightness 
profile on various parameters. 

Results. The disc typically evolves toward a "standard" steady state, where the radial surface brightness profile smoothly decreases 
with radius as r^^-^ outside the birth ring. This confirms and extends the semi-analytical study of Strubbe & Chiang (2006) and 
provides a firm basis for interpreting observed discs. Deviations from this typical profile, in the form of a sharp outer edge and a 
steeper fall-off, occur for two "extreme" cases: 1) When the birth ring is so massive that it becomes radially optically thick for the 
smallest grains. However, the required disc mass is probably too high here to be realistic. 2) When the dynamical excitation of the 
dust-producing planetesimals is so low (< e > and < i >< 0.01) that the smallest grains, which otherwise dominate the optical depth 
of the system, are preferentially depleted. This low-excitation case, although possibly not generic, cannot be ruled out by observations 
for most systems, . 

Conclusions. Our "standard" profile provides a satisfactory explanation for a large group of debris discs that show smooth outer edges 
and S B oz r"'^. Systems with sharper outer edges, barring other confining agents, could still be explained by "natural" collisional 
evolution if their dynamical excitation is very low. We show that such a dynamically-cold case provides a satisfactory fit to the specific 
HR4796A ring. 

Key words, stars: circumstellar matter - stars: individual: /? Pictoris- stars: individual: HR4796A - planetary systems: formation 



^ 1. Introduction 

.5^ 1.1. the ubiquity of ring-like features 

Dusty debris discs have been detected by their infrared ex- 
cess around ~ 15% of nearby main sequence stars (e.g. 
iBackman & Paresc3ll993l) . More than a dozen of these discs 
have also been imaged, mainly in scatte red light, since the initia l 
observation of the /3 Pictoris system bv lSmith & Terrilj (|1984|) . 
One unexpected result from these images is that almost no sys- 
tem displays a smooth extended radial profile: the usual mor- 
phology is the presence of rings (or annuli) where the bulk of 
the dust populat ion is located. This ring morphology is in fact 
so common that IStrubbe & Chiand ((2006) pointed out that the 
debris disc phenomenon could more appropriately be renamed 
debris "ring" phenomenon. Even for the archetypal debris "disc" 
yS Pictoris, which has been imaged from 5 to a few thousand AU, 
the bulk of the dust is probably concentr ated in a ra ther narrow 
region between 80 and 120 AU (e.g. Augereauetal. 2001). One 
of the few systems actually resembling an extended "smooth" 



disc might be Vega, for which Spitzer mid-infrared images show 
a rather smooth radial luminosity profile (S u et al. 200 5). At the 
other end of disc morphologies, among th e most striking ring 
features are the ones around HR4796 (e.g Javawardhana et al.l 
1998; Koerneretal. 1998; Schneider et all 1 19991). H P 139664 
dKalas et al.ll2OO60 . and Fomalhaut ( K alas et alj|2005l) . 

The characteristics that most differentiates one ring-like sys- 
tem from the other is the sharpness of the luminosity drop at the 
inner and outer edges of the rings. We shall in this paper focus on 
the outer edge issue, for which systems can be basically divided 
into two categories (see also lKalas et alj|2006l) : 

- "Sharp edge" rings, displaying abrupt surface brightness 
drops, sometimes as steep as r"^ beyond the ring. The most 
representative members of this group are HR4796A and 
Fomalhaut. 

- "Smooth edge" rings, with no sharp outer edge and a sur- 
face brightness drop beyond the ring in r^^ or r The most 
famous examples are here /3 Pic and AU Mic 
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See TablU for a list of outer-edge profiles for several resolved 
debris discs. 
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Table 1. Geometry and surface brightness profile for a selection of debris disc systems resolved in scattered light images. 



System^ 


Orientation 


Detected Radial Extent'' 


Surface Brightness S B oc r" 


Reference 


HD 53143 


~ face-on 


55 - llOAU 


a~-3 


Kalas et al. (2006) 


jSPic 


edge-on 


127 - 193AU 


a ^-3 


Golimovski et al.(2006) 






193 - 258AU 


a ^ -4 




HD 32297 


edge-on 


on 1 ATT 

80 - 170AU 


a 3.5 


Schneider et al. (1999) 






170 - 350AU 


a 3 (averaged over the 2 ansae) 




AUMic 


edge-on 


32-210AU 


a 3.8 (averaged over the 2 ansae) 


Fitzgerald et al (2007) 


HD 139664 


edge-on 


83 - 109AU 


a ~ -4.5 


Kalas et. al. (2006) 


HD 107146 


~ face-on 


130- 185AU 


a ~ -4.8 ±0.3 


Ardila et al. (2004) 


HD 181327 


face on 


100 - 200 AU 


ff ~ -5 


Schneider et al. f2006) 


Fomalhaut 


66° inclination 


140 - 160 AU 


a ~ -6.1 


Kalas et al. (2006) 


HR 4796A 


73° inclination 


70 - 120AU 


a ~ -7.5 


mhhai et al. (2005) 



" For many of these systems, many additional features, i.e., warps, clumps, etc., have been observed (it is especially true for /? Pictoris) but we 
focus here on the main issue of average radial profiles 

* The radial extents and surface brightness profiles are given for regions beyond the likely "birth ring" - the region where scattering luminosity 
peaks and where most parent bodies are believed to reside (see Sec.O 



The presence of ring features is commonly attributed 
to some sculpting mechanisms including, in partic ular, the 
presenc e of a massiv e planet (see, f or in stance, lOuillenI 
(2007a) for Fomalhaut, 'Freistet ter etal] (|2007|) for /3 Pictoris, 
[Wya ttetal. (1999) for HR4796A or the more general study 
of lMoro-Martin & Malhotral ( l2005h ). The gravitational effect of 
such a planet can truncate or create gaps in the disc, either di- 
rectly on the dust particles or indirectly on the planetesimals that 
produce these particles. 



1.2. outer edges 

Nevertheless, while planets may be a natural and easy expla- 
nation for sculpting the inner edges, outer edges are a different 
problem. They are difficult to explain in the light of one well 
established fact about debris discs, i.e. that the observed dust is 
not primordial but steadily produced by collisions, through a col- 
lisional cascade starting at much larger parent b odies, maybe in 
the planetesimal size range (e.g. Lagrange et al.ll200 0). In this re- 
spect, even if there is a sharp outer edge for the parent body pop- 
ulation, collisions would constantly produce very small grains 
that would be launched by radiation pressure onto eccentric or 
even unbound orbits, thus populating the region beyond the outer 
edge and erasing the appearance of a narrow ring over timescales 
which might be shorter than those for gravitational sculpting by a 
planet. This issue is a critical one, since these small grains dom- 
inate the total geometric cross section and thus the flu x in s cat- 
tered light (see the discussion in lThebault & Augereau H2007b . 

One possible confining mechanism for the sharp outer edge 
is gas. For discs transiting between gas-dominated phase (proto- 
planetary discs) to dus t-dominated phase ( debris discs), narrow 
dust rings may arise (iKlahr & LinI l2005l) . iBesla & Wul (l2007h 
further demonstrate that there exists an instability with which the 
residual gas collects grains of various sizes (even those subject to 
radiation pressure) into a narrow belt. However, this mechanism 
requires at least comparable amount of gas and dust. This may be 
difficult to justify for most debris discs, which are evolved sys- 
tems where the amount of gas is probably too low to prevent the 
smallest grains to be launched onto very eccentric orbits smooth- 
ing out any sharp outer edge. 

"Razor sharp" outer edges are thus very difficult to explain 
in the presence of this unavoidable outward launching of small 
grains. However, although no perfect abrupt outer edge has in- 



deed been observed Q (unlike inner edges, which are in some 
cases, like Fomalhaut, almost razor sharp), a great variety of 
outer edge profiles does exist, from relatively smooth to very 
steep (see Tab.[TJ. 

In this study, we address the issue of how these different 
profiles can be physically achieved: can this diversity be ex- 
plained by the sole "natural" evolution of a collisional active disc 
steadily producing small, radiation-pressure affected grains, or is 
(are) additional mechanism(s) needed? We consider initial con- 
ditions which are a priori the most favourable for creating sharp 
edges, by assuming a population of large parent bodies confined 
within an annulus with an abrupt cutoff at its outer edge. How 
this initial confinement may have come about is itself an in- 
teresting question but is not the focus of the current paper (see 
however the discussion in Sec. 15. 2b . The outcome we consider as 
a reference for our investigation is the radial surface brightness 
(SB) profile in scattered light, since this is an observable which is 
reasonably well constrained for most imaged debris discs (either 
directly observed or obtained by de-projection). We consider the 
nominal case of a disc seen edge-on, but results can easily be 
extrapolated to face-on systems, since SB profiles for both ori- 
entations tend towards the same radial dependence far from the 
birth ring (given the same radial dust distribution). 

In several previous studie s , all a ddressing the s pecific 
p Pic case (Lecavelier etakl 119961 : lAugereau et alj 120011 : 
iThebault & Augereau I l2005h . it has been argued that the "nat- 
ural" SB profile outside the collisionaly active parent bodies 
belt, or "birth ring", falls off as ~ r^^. This is based on the 
assumption that all particles produced in the birth ring have a 
size distribution which scales as dN/ds oc s^^-^, (as expected 
for an ide alized infinite collisional cascade at equilibrium, see 
iDohnanvT 1969), down to the radiation blow-out limit s - 
(where the ratio of radiation pressure to gravity /3 = 0.5). The 
smallest radiation pressure-affected grains, which dominate the 
light receiving area, are then diluted along their eccentric or- 
bit. This geometri cal spread results in SB oc r^^. However, 
IStrubbe & Chiand (2O06) argued that, since high-/? grains spend 
a long time in the collisionaly inactive region beyond the birth 
ring, the dN/ds oc s^^-^ collisional equilibrium law should only 
apply to the small fraction / of these grains which are present 



' HD141569A could be a possible exception, but this system is 
likely not a true debr is disc, being significantly younger and gas rich 
iJonkheid et alj|2005) . but a member of the loosely defined "transition 
object" category. 



Thebault and Wu: Outer Edges of debris discs 



3 



in the collisionaly active birth ring. This results in a 1// excess 
of the disc-integrated number of small grains, which in turn re- 
sults in a flatter S'Bprofile in r"^^. They applied their theory 
to the AU Mic disc and reproduced the obs erved SB profile, 
spectr al energy distribution and disc colour .Strubbe & Chiand 
(120061) further argued that the observed SB profile depends only 
weakly on the radial and size distributions of grains within the 
birth ring. The discs which exhibit a fall-off" sharper than r"^^ 

are thus puzzling in the face of this theory. 

The innovative model of IStrubbe & ChiangI (l2006h is build 
on analytical derivations and Monte-Carlo modeling which did 
not actually treat the collisional evolution of the system and re- 
lies on several simplifying assumptions. The main one is that the 
size distribution is fixed and is assumed to follow the idealized 
dN/ds oc s^^'^ scaling (corrected by the fraction 1//), whereas 
several st udies have shown that this la w cannot hold in real sys- 
tems (see lThebault & Augereaull2007 ,and references therein)0. 
Another issue is that when evaluating collisional life-times, only 
the vertical velocity of the grains was taken into account, thus ne- 
glecting their radial movement which can be appreciable, if not 
dominant for the smallest grains. Finally, the specific dynamics 
of the small radiation-pressure-aff'ected grains, in particular the 
fact that they suffer much more frequent collisions and at much 
higher velocities, is not taken into account. 



depend on average values of the encounter velocities Av^iajja'j')- 
For objects not subject to significant radiation pressure, relative 
velocities are obtained from the eccentricity and inclination dis- 
tributions t hrough the classical expre ssion, valid for randomized 
orbits (e.g. lLissauer & Stewarll993h : 



1/2 



(1) 



where Vkep{ia) is the Keplerian velocity at radial distance rj^. 

For the smallest particles, the effect of radiation pressure, 
which places objects on highly eccentric orbits, is taken into ac- 
count. Inter annuli interactions, induced by the significant radial 
excursion of these bodies are considered, and Av(,fl,,;,a',,v) are de- 
rived through separate deterministic N-body runs. 

Collision outcomes are divided into 2 types, cratering and 
fragmentation, depending on the ratio between the specific im- 
pacting kinetic energy E^oi and the specific shaterring energy 
2*, which depends on object sizes and composition. In both 
regimes, the size distributions of the newly produced frag- 
ments are derived through detailed energy scaling pres crip- 
tions, which are presented at length in Thebault et al.l (l2003h and 
iThebault & Augereau I (12007 ). Possible reaccumulation onto the 
impacting objects is also accounted for 



2. Our Approach 

We re-address these issues using a numerical approach quantita- 
tively following the collisional evolution of the full system. We 
start with a birth ring of parent bodies in a perfectly confined 
annulus and let it collisionally evolve. The temporal as well as 
spatial evolution of the size distribution are followed, taking into 
account the radial excursions of high-y6 particles. As previously 
mentioned, we derive for each simulation the surface brightness 
profile in scattered light. 

2. 1 . numerical model 

We use a statistical particle-in-a-box model to follow the evolu- 
tion in size and spatial distribution of a population of collision- 
ally interacting bodies. This code has initially been developed, in 
its single-annu l us ver sion, for the study of the inner /3 Pic disc 
(IThebault et al.l l2003l) . and later upgraded to a multi-annulus 
version (i.e., with 1-D radial resolution) for the study of colli- 
sional processes in extended debris discs (Thebault & Augereau] 
l2007h . The detailed description of the code can be found in these 
two papers, and here we recall some of its main characteristics. 

The entire system is spatially divided into concentric ra- 
dial annuli. Within each annulus of index ia, the solid body pop- 
ulation is divided by size into bins that cover a broad size 
range spanning from kilometre to micron. With a standard log{2) 
size increment, this requires ~ 100. Evolution of the particle 
number within one (ia, i) bin (/ being the size distribution index) 
is contributed by all destructive impacts between (ia, i) objects 
and bodies from other (ia', /') bins as well as all impacts between 
other size bins producing new (ia, i) objects. Collision rates are 
estimated statistically. These rates, as well as collision outcomes, 

^ AU Mic is an M star with weak radiation but where stellar wind is 
believed to act on small grains in an equivalent way as radiation pressure 
does around more massive A stars (e.g lAugereau & Beustii2006l) . 

^ The need for realistic size distributions departing from fixed pow er 
laws has been very recently emphasized by iFitzgerald et al.l dlOOTl) in 
their analysis of the AUMic data 



2.2. Set Up 

Our numerical model requires the following inputs: the ring's av- 
erage distance from the star, r^R, its radial width Ar^R, the aver- 
age free eccentricities (e) and inclinations (/) of the parent bodies 
(non affected by radiation pressure) and the initial particles' size 
distribution and total mass. We chose to parameterize the latter 
by Mjiis,, the total mass of objects with sizes s < = I cm, even 
though our simulations include bodies up to 10 km, because the 
total "dust" mass is a parameter which can often be constrained 
from observations. We also assume equipartition between eccen- 
tricity and inclination, so that (e) - 2{i), at least for the parent 
bodies. 

Due to the numerical cost of the detailed size distribution 
evolution procedure, (e) is assumed fixed, independent of posi- 
tion and time. The first independency is justified as we are con- 
sidering a relatively narrow birth ring. The time independency 
is an acceptable simplification when considering that the eccen- 
tricity disper sion is imposed by the largest bodies present in the 
system (e.g. IOuillenll2007bt) and that, according to our cuiTent 
understanding of debris discs, i.e. systems in which the bulk of 
the planetesimal accretion process is already over, these bodies 
should be large planetary embryos that are too big and isolated to 
be significantly afifected by collisional erosio n over the ~ 10^ - 
lO^yrs timescale considered here (see, e.g. IWvatt et al] 120071; 
iLohne et al.ll2007l) . This is also why these dynamics-imposing 
embryos are left out of the collisional cas cade numerically stud- 
ied h ere (for more on the subject, see IThebault & Augereau I 
' 2007[ as well as the discussion in Sec l5.1.2l l. 

For the sake of clarity, we consider a nominal case (Tab|2]i 
with set-up matching as closely as possible the /3 Pictoris system, 
i.e. Mdusi = O.lMe and (e) = 0.1 (e.g. Augereau et al .11200 lb . 
For the radiation pressure blow out size, we take so.5 = 5//m,the 
value derived for compact silicates around a j3 Pic like A5V star 
The system is divided into two radial zones: 

- The parent body zone, or "birth ring", located between 
and Tout with the center at tbr and divided into 6 annuli, 
where we follow the collisional evolution of the whole solid 
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Table 2. Nominal case set up. The fields marked by a * are ex- 
plored as free parameters in the simulations. 



Radial extent of birth ring* 


80 < r < 120 AU 


Number and radial width of annuli* 


6X6.66AU 


Initial surface density within ring 


uniform 


Total "dust" mass M^^^t < 1 cm)* 


0.1 Me 


Stellar spectral type* 


A5V 


Radiation Blow-out size* 


io.5 = 5 yum 


Size range modelled 


2fim< s <10km 


Number of size bins 


99 


Initial size distribution' 


dN oc s-^ ^ds 


Dynamical excitation* 


<e>o =0.1 =2</>o 



body population, from a maximum size imax = 10 km that 
sits in the planetesimal size range to a minimum size s^in - 
2jjm below the radiation blow-out limit. The initial size dis- 
tribution in the entire range is assumed to follow the ideal- 
ized coUisional " equilibrium" distribution, dN(s)o oc s^^'^ds 
( lDohnanvilll969 l). Our runs show that this choice is not cru- 
cial: in the relevant dust-size range, the size distribution is 
quickly relaxed toward a new steady state with a profile that 
is independent of the initial choice and that deviates signif- 
icantly from a Dohnanyi-like power law (see section |4] and 
Fig©. 

The outer zone, which is devoid of particles at the beginning 
of the runs and gets progressively populated by small grains 
coming from the birth ring. Consequently, we only follow 
here grains with s < s^^^, where s"^^^ roughly corresponds 
to the biggest grains able to leave the parent body region 
and is taken conservatively to have a radiation-to-gravity- 



ratio of /3(s', 



0.1. Spatially, this region is divided into 



3 annuli just outside the main ring plus one additional, in- 
finitely extended "buffer" annulus. Within the latter zone, no 
collisional evolution is modeled, and only the orbital evolu- 
tion of the grains is considered: either escape of the system 
for unbound grains or progression to the apoastron and re- 
turn to the inner annuli for the bound ones. The radial ex- 
tent of the outer zone (not including the "buffer" annulus) is 
set to ~ 4-rBR, which is typically the extent of the"outer" re- 
gion considered for the 2 most famou s birth-ring/outer-zone 
systems, i.e., 13 Pic and AU Mic (e.g. lAugereau et al.ll200ll 
IStrubbe& Chiang|l2006h . 

The radial surface brightness profile in scattered light is then 
synthetically computed using the dust size and radial distribu- 
tions, assuming an r^^ dilution of the stellar flux and isotropic 
scattering (although different scattering properties are also ex- 
plored). 

We explore around the nominal set-up (see Tab|2]i, especially 
for the two fundamental parameters which are M^ust and the dy- 
namical excitation of the system (as parameterized by (e)). 



3. Numerical Results 

3. 1 . nominal case: smooth edge and r^-^ profile 

FiglUpresents results obtained for our /3 Pic-like nominal case. 
As can be clearly seen, the system rapidly (~ lO^yrs) reaches a 
steady state after which the SB profile no longer evolves signifi- 
cantly. We stop the integration at lO^yrs. This ~ lO^yrs timescale 
is the time it takes for the small grains that fill the outer radial 
zone to reach an equilibrium between collisional production and 
destruction. A direct consequence of this fast evolution is that 
the initial sharp outer edge of the parent body region is quickly 
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Fig. 1. Evolution of the midplane luminosity profile in scattered 
light. Nominal case: parent body ring [80, 120]AU, Mdust = 
0.1 Me, (e) = 0.1. Each midplane luminosity has been renormal- 
ized to its value at 50 AU. The dark grey area represents the par- 
ent body region where all mass is initially located. The narrow 
light grey area represents a r"^^ slope reconstructed backwards 
from the final luminosity value at 400 AU with a +15% width. 
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Fig. 2. Nominal case at f = lO^yrs. Respective contributions of 
different grain populations (parameterized by their /3 value) to 
the total scattered flux, as functions of the radial distance. 



smoothed out. Once the steady state is reached, the profile in the 
r > rout region lies very close to S B(r) oc r"^^. 

When looking at the respective contributions from different 
dust populations to the total profile, it appears that the scattered 
flux is, in the 1 20^00 AU region, dominated by high-/? grains 
in the 0.25 < /3 < 0.4 range (Fig|2]i. These grains have or- 
bital eccentricities in the 0.33-0.75 range and apoastron in the 
240-700AU region and thus spend a large fraction of their orbits 
in the domain located between rout and 4rBR. Grains with even 
higher /? (close to 0.5) only weakly contribute to the flux in the 
'"out to 4rBR region because they have orbits whose apoastron is 
often at several lOOOAU and will thus spend most of their time 
outside the < 400 AU region considered here. 
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Fig. 3. Same as Figd] but with Mdust = O.OOIM®, all other pa- 
rameters being the same. Note that the final timescale is here 
t - lO^rs. 

Table 3. Results obtained for different values of the system's 
dynamical excitation (e) {-2 (/) ) and width of the birth ring 
ArBR B Also shown are results for the nominal run but with the 
SB computed using a different scattering phase function (see text 
for details) 



Run Outer edge sharpness^ SB profile^ 







5B(ro„t)/5B(ro„t+io) 


slope a 


nominal case 




2.50 


-3.51 


<e>= 0.2 




2.44 


-3.48 


<e>= 0.035 




2.78 


-3.52 


<e>= 0.01 




3.33 


-3.49 


ArsR = 20AU 




2.05 


-3.51 


ArsR = lOAU 




1.96 


-3.54 


scattering anisotropy ; 
scattering anisotropy j 


1=0.5 
i=0.8 


2.22 
2.22 


-3.68 
-3.98 



" when varying its outer edge Vant, its inner edge being fixed at 80AU 
* luminosity drop at the outer edge of the ring, as measured by the 
flux ratio between the outer edge of the birth-ring and lOAU beyond it 
average value in the i.SrsR- 4rBR region 



For a case with a much less massive {Mjust - 0.001 Me) disc, 
we see that the system reaches a similar outer SB profile (SB 
oc r"^-^ outside rout) but with a much longer timescale: a few 
10^ yrs instead of ~ 10^ yrs in the nominal run (Fig.[3]i. Similar 
SB profile is also obtained for a higher mass run (Mdust - IMg). 
Once again, the only difference with the nominal case is the pace 
at which the steady-state is reached (almost 10'*yrs here). 

We further explore the parameter dependences by running 
a series of simulations varying the width of the birth ring and 
the dynamical excitation of the system. We also investigate the 
importance of the scattering function assumption by explor- 
ing anisotropic cases, assuming a iHenvev & GreensteinI (1 194 lb 
phase function and changing the asymmetry factor | g \. The re- 
sults for all these runs are summarized in Tab|3] showing the 2 
main outcomes of interest for the present problem: the luminos- 
ity drop at the outer edge of the ring and the average slope of 
the luminosity profile in the l.SrBR^rBR region. The nominal 
results are robust: the luminosity drop at the ring's outer edge 



Fig. 4. Same as FigH] but with a very high dust mass of Mdust - 
lOMe. 



is always comprised between 2 and 3, whereas for the slope we 
get a - -3.51 + 0.03. The only way to reach a somehow steeper 
slope is to assume a strong anisotropy of the scattering function. 
But even for a rather extreme g=0.8 case, a ~ -3.98. 

3.2. "extreme cases" producing different profiles 
3.2.1. very massive disc 

In order to push our parameter exploration to the limits, we keep 
increasing the total initial dust mass Mdust until we observe a 
departure from the nominal case SB profile. Such a departure is 
observed for Mdust - lOM® (FigHJi. This extremely massive case 
distinguishes itself from the others by the fact that a sharp lumi- 
nosity drop is maintained at the outer edge for ~ lO^'yrs. What 
keeps the outer regions at a very low luminosity is the fact that 
the parent body ring is so densely populated that grains pushed 
outward by radiation pressure cannot freely escape it without 
experiencing a collision. This can be illustrated by looking at 
a simplified parameter, i.e. the geometrical radial optical depth, 
defined (at a given distance from the star ro) by 

T„d(ro)=2. ^ ^2) 

ia=w 

where dNia/ds(r) is the differential number of i-sized particles 
in a radial annulus of index ia centered at radial distance r, H 
is the vertical height of the disc at that distance and /O is the 
annulus index corresponding to the radial distance ro. Trad is of 
course simply the optical thickness of the disc to stellar photons. 
Since particle orbits are never straight radial lines, this quantity 
is only a first approximation of their real in-plane (or horizontal) 
collisional probability. However, for the smallest grains this is a 
relatively good first-order approximation (see the more thorough 
discussion on horizontal and vertical collision probabilities in 
Sec.lSJil. 

Fig|5]shows that Trad ^ 1 over most of the source rin^in the 
early epoch, so that few grains can escape the birth ring without 



When this is the case, one should also include the attenuation of 
stellar light over the disc to be self-consistent. However, we believe 
doing so would not qualitatively alter our conclusion. 
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Fig. 5. Evolution with time of the radial optical depth (defined as 
in eq. ||2l) for the very high dust-mass case Mdust = IQM^. 
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Fig. 6. Same as FiglT] but for a dynamically "very cold" system 
with (e) = 0.001. 



colliding with another grain, hence the density (and luminosity) 
depletion in the r > regions. However, this radial optical 
depth steadily decreases over time, due to rapid mass erosion by 
energetic collisions within the birth ring. It drops below unity af- 
ter ~ 10^ yrs, and is ~ 0.6 by the time (~ 10^ yrs) the sharp outer 
edge is smoothed out (see FigH]). After that, the system behaves 
like the nominal case. We re-examine this radially optically thick 
case in more detail in Sec l5.ll 



3.2.2. dynamically cold system 

Another way of obtaining a departure from the standard SB pro- 
file is to decrease the dynamical excitation of the system to a 
very low value, typically (e) < 0.01. Fig|6] shows the results 
for the case (e) = 0.001 ( - 2(/)), where a sharp outer edge is 
maintained outside the birth ring for the entire duration of the 
simulation (10^ yrs). 

The abrupt luminosity fall-off at the outer edge is here due to 
a global depletion of small grains directly resulting from an im- 



Fig. 7. Geometrical cross section per logarithmic size range, in- 
tegrated only for grains inside the birth ring, for systems with 
different dynamical excitation. These curves at obtained at the 
end of the integration (10^ yrs) and are normalized to have the 
same value at the large grain size. 



balance between the collisional production and destruction rates 
of small and large particles in such a dynamically cold disc. This 
can be understood in the following way: 

For a given dust mass, a lower dynamical excitation does 
not change the collision rate between large objects not aff'ected 
by radiation pressure. Indeed, because of the equipartition (e) = 
2(/) the decrease in (Av) is exactly compensated by the increase 
of the particle number density due to the reduced thickness. 
However, the lower (Av) values mean that collisions will be less 
destructive and produce less smaller fragments. This means that 
the rate at which small grains (the ones significantly affected by 
radiation pressure) are produced is significantly reduced. On the 
contrary, the rate at which these small grains are destroyed is 
higher. Indeed, collisions velocities for impacts involving them 
are not significantly reduced by the small (e) for parent bodies, 
since small grain dynamics is predominantly imposed by radia- 
tion pressure. Furthermore, the rate at which such impacts occur 
is increased, compared to the nominal case, because of the in- 
creased radial optical depth (see Sec l5.1.2l for a more detailed 
analysis). 

As (e) decreases, this imbalance becomes more severe and 
the depletion of the smallest grains is more acute. This is illus- 
trated in Fig|7]displaying, for different values of the system's dy- 
namical excitation, the respective contributions of different grain 
sizes to the total geometrical cross section cr. For the nominal 
case ((e) = 0.1), we obtain the standard result that cr is dom- 
inated by the smallest grains close to the cut-off size. As (e) 
gets smaller, however, the contribution of these smaller grains 
progressively decreases. Below the limiting value (e) ~ 0.01, 
this effect is so pronounced that the system's optical depth is 
no longer dominated by the smallest grains, but by much bigger 
particles in the lOO-lOOOio.s range. These particles have their or- 
bits largely confined within the birth ring. This explains why a 
sharp luminosity decrease is observed at the outer edge of the 
birth ring. 

Note that contrary to the very high mass case, the sharp outer 
edge does not smooth out with time but persist throughout the 
lO^yrs of the simulation. 
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4. Analytical derivation: The 'universal' 



-3.5 



profile 



Our numerical exploration has shown that the r^^-^ surface 
brightness profile beyond the outer edge seems to be the most 
generic outcome for a collisional ring system under the action 
of stellar radiation pressure. We reinvestigate this issue from an 
analytical point of view and derive simplified formulae confirm- 
in g this result. We take her e as a basis the analytical approach 
oflStrubbe & ChiangI (l2Q06|) and extend it to more general cases 
regarding grain size and spatial distributions. 

We firstly assume that, within the birth ring, particles follow 
a power-law size distribution of index q (instead of fixing q - 
-3.5 for a Dohnanyi equilibrium): 



c/A^BR °^ s''ds. 



(3) 



Since these grains spend most of their orbits in the empty re- 
gion outside the birth ring, their total number integrated over fhe 
whole system (Mot) will be boosted by a factor l//(e), where 
/(e) is the fraction of an orbital period a b ody of eccentricity 
e ~ B /(l - /3) spends within the birth ring (IStrubbe & Chian j 
120061) . so that 



fie) /(e) 



(4) 



We also followlSttubbe & Chiang' ('2006') in making the sim- 
plifying but reasonable assumption that all high-/? grains are on 
average mostly seen near their apoastron. Thus, at each given 
distance r from the star, the optical depth is dominated by parti- 
cles of size 



Sdomir) 



1 



1 



iO.5- 



(5) 



Since we are interested in the region r < 4rBR, only grains with 
their apoastron 



fl(l +e) = -^-^11 + 



l-2j0\ l-p 



'"BR = 



^BR 

1-2/? 



<4rR 



(6) 



are important. This corresponds to yS < 0.4 and thus e < 0.67. 
Note that, in their analytical derivation (this assumption was 
relaxed in their Monte-Carlo model), Strubbe & Chiang (2006) 
only considered extremely eccentric e ~ \ grains, which spend 
most of their orbits outside the considered < 4rBR region. 

To estimate the enhancement ratio 1 //(e), we do not adopt 
the asymptotic expression (valid for eliml) /(e) oc (1 - e)'-^ 
taken by Strubbe & Chiang (2006) but instead derive this value 
from Kepler's equation for a typical particle produced at the mid- 
dle of the birth ring: 



/(e) oc (£2 - £1) - e(sin(E2) - sin(Ei)) 
where 



£, = 



and 



E2 — acos 



i ll _ ^BR + ArBR 



(7) 



(8) 



The vertical geometrical optical depth then scales with radius as 



1 dN,ot(Sdom) 2 / X 

^^^'^ " ; dr '^'""^''^ 

1 1 dNBRiSdom) 2 



,('•)■ 



/(e) r dr 

Using the radial dependence of Sdom given by eq|5] we obtain 



(9) 



dr 



(, rBR\ 'I ^ rsR 

V-—) — 



(10) 




dNBB=s Ms 
assuming dNt^t=s"^^d 



distance to tlie Star (in Fhk units) 

Fig. 8. Radial profile of the vertical optical depth, obtained with 
our analytical fit of eq[TT] assuming q = -3.5 (with a ring width 
ArBR - 0.1 tbr). As a comparison, we also plot the profile ob- 
tained with the (incorrect) assumption that the size distribution 

is dNtot s^^'^ds in the whole system. 



Inserting Eqnsl4ll5]and[T0linto Eq|9]leads to 



Tj.(r) oc 



-q-A 



fBR 



(11) 



where /(e) is given by solving Equ|2] 

As a test, we first consider the same fiducial case as 
IStrubbe & Chiang (2006), i.e. the Dohnanyi q- -3.5 value. The 
radial profile of Ti(r) obtained this way scales as r"'^ (Fig. [8). 
Departures from this slope are relatively limited, even in the 1- 
1.5 rBR region, where Tj.(r) has a ~ ^175 dependence. If we 
apply the usual rule of thumb that for an optical depth profile 
T^(r) oc r", the midplane SB scales as r''"^(e.g. Nakano 19900, 
we obtain S B(r ) oc r^^-^ in most of the o uter region, confirm- 
ing the result of lStrubbe & Chiang! (l2006h . For comparison, we 
also plot on the graph the T±(r) profile derived when (incor- 
rectly) assuming that the Dohnanyi law ap plies to the entire sys- 
tem, i.e., dNtor oc s^^^^ds (as was done in .Lecavelier et al.lll996t 
lAugereau et al.l2001l:lTTiebault & Augereau Il2005h . Not surpris- 
ingly, we recover the asymptotic dependence T±{r) oc r"^, corre- 
sponding to the S B(r) oc r^^ slope derived in thes e past studies. 

However, as has been often pointed out (e.g. iThebault et al] 
2003), size distributions follow a power law with index q = -3.5 
only in unrealistically idealized systems, with an infinite size 
range and no size dependence for collisional processes. In our 
simulations, the steady state size distribution within the birth 
ring significantly departs from this value (see Fig|9]l. This de- 
parture from the Dohnanyi law, especially in the crucial size do- 
main of the smallest grains, is a well know results which has 
been obtained and discussed in several previous studie s (e.g. 
(Tampo Bagatin et al. 1994; Thebauhetal. 2003; Krivo v et all 
HooilThebault & Augereau l2007l:lLohne et al.l2007 ). The main 
reason for this behaviour is that radiation pressure introduces 
a natural minimum cut-ofl" sq.s in the size distribution, so that 
bodies of size si just above io.5 are overabundant because of 
the lack of small s < io.5 impactors that could destroy them. 
Subsequently, larger bodies, of size S2, which can be destroyed 



' This assumes grey scattering and tiiat the disc's vertical thickness 
scales as H oc r. 
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size (in S05 units) 



Fig. 9. Grain size distribution inside the birth ring (dNBR/ds), for 
all cases presented in Sec|3] The size range here covers that of 
grains dominating the optical depth in the outer region. In order 
to facilitate comparison, all curves have been renormalized to 
their peak value at s = l.2s(i5. The thick grey line shows the 
theoretical Dohnanyi profile cINbr s'^-^ds. 



dK„=s"*ds 




10 



distance to ttic Star (in Fqi, units) 



Fig. 10. Same as Fig|8]but for different size distributions index 
q. All size distributions lead to Tj_ oc r"'-^ at large distances, but 
steeper size distributions approaching this asymptote earlier 



by impacts with si ones, are underabundant, which in turn leads 
to an overabundance of objects of size 53 > S2, etc. The resulting 
size distribution displays a pronounced wavy -pattern observable 
also in Fig|9l For the size range which is here of special interest, 
i.e. grains in the 0.15 < /3 < 0.4 (1.25^0.5 < ^ < S.Sio.s) range, 
the distribution is steeper than a Dohnanyi one and is close to 
q ~ -4. Interestingly, putting this value into Eq[TT] leads to a 
profile that is again well (in fact, better) approximated by 
the r slope (Fig.fTOll. This is because, in eq[TTl the r depen- 
dence arising from the size distribution term (1 - r^R/ry^* is 
swarmed by that from the geometric term (/"br/ (f{e)xr^) at large 
r. Only the rapidity at which the r"'^ asymptotic behaviour is 
reached depends on the index q (Fig.fTOl). 




1 - 
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Fig. 11. Ratios between the horizontal and vertical optical depth 
plotted as a function of grain size for different values of dynami- 
cal excitation. The straight lines associated with each case mark 
the values of Tiad/Tj. (see eq. |2]i. The grains are assumed to be 
launched with local Keplerian circular velocity at the centre of 
the ring. The fractional ring width ArsR/r^R takes the nominal 
value of 0.3. Decreasing this value will move the peaks towards 
the right. 



5. Discussion 

5. 1 . How to escape the universal r^-^ profile 

Our numerical explorations show that there are only 2 ways 
not to end up with the standard result with no sharp edge and 
S B oc r"^^: a very massive, radially optically thick disc or a dy- 
namically cold system. For these 2 "extreme" cases, the analyti- 
cal derivation of the previous section becomes invalid. Although 
the reasons why the analytical study no longer holds are different 
for each case, they are nevertheless for both cases distinct con- 
sequences of the same crucial characteristics of the system's dy- 
namics when pushed to its limits: the strong imbalance between 
the collision rates and velocities of large and small particles. 

Therefore, before discussing these 2 cases in more detail, 
let us present the mechanisms at play behind this large/small 
particle dichotomy. 

For large bodies, we have seen that collision velocities can 
be derived using the standard expression for their dynamical rms 
excitation (Eq.l). Since for these objects equipartition between 
in-plane and off -plane motions results in a fixed (e)/(/) ratio, it 
follows that collision velocities are directly proportional to the 
system's inclination whereas collision rates are independent of 
(/). That is to say that the collisional optical depth is directly 
proportional to the vertical optical depth Tj_ (which does not 
vary with (/)). This is no longer the case for small grains pushed 
on eccentric orbits by radiation pressure, which can sample the 
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whole radial extent of the birth ring. Radiation pressure induced 
eccentricities are equal to 



= 1- 



(l-2fl/?/r)(l-(g^)) 
(1 



1/2 



(12) 



where a is the semi-major axis of the parent body producing the 
small grain, and r the radial location of its release. For these ob- 
jects, in-plane motions become much more important than verti- 
cal ones. This can be quantitatively estimated by computing the 
in-plane collisional optical depth r//. This quantity is obtained 
by integrating, for a given particle, over distances that are trav- 
elled by this grain relative to its background. 



Vo|t/f no". 



(13) 



where v is the grain velocity, vo the mean velocity of particles in 
the region, which we approximate to the local circular Keplerian 
velocity, n and cr are the local number density and grain cross 
section, and their product is related to the vertical optical depth 
by Tj. — 2Hncr. The integration {j>) is over the whole orbit and 
n<T is non-zero only inside the birth ring. For a given particle 
T// depends on 2 main parameters: its orbital eccentricity and, 
for the smallest grains, the birth ring's radial extent (since this 
is the only region where collisions can occur). For large objects, 
equipartition between {e) and (/) means that r// is simply equal 
to (see FigfTTTl. which is another way to say that impacts in the 
in-plane and vertical directions are equally important. As grains 
get smaller, however, the ratio t///tj. increases, reaching a max- 
imum when the grain's orbit has the same radial extent as the 
birth ring, hence the peak at s ~ 7 so.5 for the nominal birth ring 
considered in FigfTTI 

The imbalance between small and large grains collision rates 
is directly proportional to epl{e). Since ep only weakly depends 
on (e) (EqfTSli. this means that the imbalance increases for dy- 
namically cold systems, as clearly appears in FigfTT] 

5.1 .1 . radially optically thick disc 



In Sec l3.2.ll we identified one consequence of the boosted col- 
lision rates of small particles, which is that, for very high Mdusi 
values, high-yS particles have a significant chance to experience 
a collision before escaping the birth ring. In this situation, the 
analytical derivation of the previous section is invalid because it 
implicitly assumes that all high-/? grains produced in the birth 
ring eventually reach their apoastron far in the outer regions. 

To characterize this collisional optically thick case, we used 
in Sec |3.2. li the parameter t^^^, as defined in Eq|2] and empiri- 
cally found that the limiting Trad value is around 0.6. The validity 
of the Trad parameterization is confirmed by FigHT] which shows 
that this simplified parameter is actually a good approximation 
of the "real" in-plane (or horizontal) collisional optical depth r// 
for small grains close to the .so.5 limit. We thus keep this easily 
defined parameter (as compared to t//) as a good approximate 
limiting criterion. Writing as a first approximation that 



Trad 



Ar, 



BR 



1 



'"BR <e) 



<Tx> 



BR, 



(14) 



we see that there are two ways of reaching high Trad values: ei- 
ther by increasing the total dust mass of the system (as parame- 
terized by Tx or Mdiisi) or by decreasing its dynamical excitation. 

We have seen in our numerical exploration that only playing 
on the Mdusi (or t^) parameter requires to reach values of the 



order of lOM® for our nominal case of a disc with (e) = 0.1 
and ArBR/rBR - 0.3. However, our simulations show that global 
collisional erosion in such a high-mass system is so intense that 
the Tiad ^ 0.6 regime can only be maintained for ~ lO^yrs be- 
fore the birth ring becomes radially transparent for high-/? grains, 
at which point the sharp outer edge vanishes (see FiglDl. One 
might argue that the largest bodies in our simulation are in the 
10 km range and thus can be quickly eroded in this high density 
case, while realistic systems may have a large mass reservoir 
contained in larger bodies and could sustain a longer period of 
extreme dustiness. However, in this case the total mass of solids 
in the system, which is already M,ot ~ 10^ Me in our high-mass 
run with s^ax = 10km, would reach unrealistically high values. 
Another strong argument for ruling out this high-mass case is 
that the presence of lOM® of dust is not supported by observa- 
tions. Indeed, no debris disc around a main sequenc e star seems 
to contain more than ~ IMe of dust (Greaves 2005), or has a 
vertical optical depth exceeding ~ 10"^ ( Kalas 2005 ). 

The other alternative for reaching high Trad values, i.e. re- 
ducing the dynamical excitation, would require, for a typical 
Mdusi - 0.1 Me disc, (e) values as low as ~ 0.001. However, 
for such low values of the excitation, we have seen that the sys- 
tem's collisional evolution and size distribution is already in the 
qualitatively different "dynamically cold" regime, where small 
particles cease to dominate the optical depth (see next section). 

5.1 .2. dynamically cold disc: depletion of small grains 

This dynamically cold mode, with typically (e) < 0.01, is in fact 
the second way to maintain a sharp outer edge (see Sec l3.2.2t . 
In this case, the imbalance between small and large grains' col- 
lisional behaviour gets very pronounced (as it appears clearly in 
FigfTTTl. up to a point where it causes a depletion of small high- 
/3 grains. We have seen in Sec l3.2.2] that this is because small 
grains are predominantly created by impacts involving large ob- 
jects, which get less efficient for low (e), but predominantly 
destroyed by collisions involving themselves, whose efficiency 
(imposed by radiation pressure-induced motions) only weakly 
varies with (e) of the parent bodies and whose rate strongly in- 
creases for decreasing (e). As a consequence, contrary to the 
nominal case the system's optical depth and luminosity are dom- 
inated by larger ~ lOO^o.s grains that do not leave the birth ring, 
hence the sharp, and long-lived drop at its outer edge0. 

However, we are here confronted with two issues. 

The first one regards the total brightness of such discs. 
Indeed, a disc deprived from its smallest grain population (in 
the < lOOio.5 ~ 0.5mm range) should appear much dimmer 
than a disc in which these grains dominate the light scattering 
area. From Fig|2] we see that, for the same total mass of dust 
(always contained in the biggest grains), a small-grain-poor sys- 
tem with (e) < 0.01 is at least 10 times less luminous than the 
nominal case. It follows that, to reach the same fractional lu- 
minosity (assuming it is proportional to the optical depth), the 
dynamically cold system has to be at least 10 times more mas- 
sive. This has no direct implications on observed dust mass es- 
timates, which are usually derived from sub-mm or millimetre 
observations in thermal emission, but should affect the corre- 
lation between Md„sr and total fractional luminosities fi. In this 
respect, dynamically colder (or thinner) discs should have higher 
Mdusi/ fii^i^os. However, the uncertainties affecting both param- 
eters' estimates are probably high enough to accommodate a fac- 



* The analytical derivation of Sec.|4]is invalid here, since it implicitly 
assumes that small high-/? grains dominate the total optical depth 
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tor 10 uncertainty in their ratio (see for example the attempt at 
connecting these two quantifies for the specific case of AU Mic 
performed by Augereau & Beust 2006), so that our result is here 
probably not very constraining. 

The second and probably crucial issue is how likely it is to 
find such dynamic ally cold discs, with (e ) and (/) lower than 
^ 0.01. As noted by Thebault & Augereau (2007): "the only ob- 
servational constraint [on the disc's dynamical excitation] comes 
from measuring the disc's vertical thickness and deriving esti- 
mates of orbital inclinations, but such constraints are scarce". 
Edge-on discs represent the most favourable cases since H/r can 
be directly measured. However, even for the two most studied 
discs, only partial information is available. Krist et al. (2005) 
find Hjr < 0.04 for AUMic (and Hjr < 0.02 close to the po- 
sition of maximum surface density). The Pictoris disc appears 
thicker with H/r ratios as large as ^ 0.1 (Golimowski et al. 2006). 
Q. The modelling and inversion of scattered light brightness pro- 
files of inclined, ring-shaped discs do not provide many more 
constraints. The HD181327 rings for example, might have H/r 
ratio as large as about 0.1 at the positions of maximum sur- 
face density, but the actual ratios could be two times smaller 
(Schneider et al. 2006). 

There are thus large uncertainties, but it seems however that 
the 0.01 to 0. 1 range is the most realistic one for (/) (and thus {e) 
if assuming equipartition). This 0.01-0.1 range does also make 
sense when considering simple theoretical arguments regarding 
the sizes of the biggest objects within the disc. Consider indeed 
a belt of planetesimals sitting at ~ 100 AU from the central star 
If the maximum size is limited to, say, ~ 10km, then excitation 
by mutual viscous stirring leads to values of the order of their 
surface escape velocity (~ lOm/s), and collisions among them 
generate particles with dispersion velocity that does not exceed 
the same velocity. The resulting HjR ~ Vesc/v'kep ~ 0.003 and 
is dynamically cold. On the other hand, if the maximum size is 
~ 500km, then H/R ~ 0.15 and the resulting debris disc is dy- 
namically hot. This dynamically hot case probably makes more 
sense within the frame of the standard planet formation scenario, 
in which debris discs correspond to systems in which the bulk of 
planetesimal accretion proces s is already over and large plane- 
tary embryos are present (e.g. lKenvon & Bromley |[2005h . 

5.2. Application to real systems 

As a consequence, we expect our nominal result, valid in the 
{e} > 0.01 range, to correspond to the "natural" collisional 
evolution of most discs when left to themselves. This nominal 
model, with no sharp edge and luminosity falling as r^^-^ com- 
pares well against the outer region SB profiles (beyond ^ 40AU 
and ^ 120AU respectively, see TablefTli for the two perhaps most 
emblematic debris discs: /5 Pic and AU Mic. For these two sys- 
tems, which are amongst the few ones for which (partial) disc 
thickness estimates are available, this result is in good agreement 
with the observed (/), which is for both discs in the > 0.01 range, 
and thus in the dynamically "hot" regime displayed in FigfTFI 
For AU Mic, this conclusion confirms the initial results obtained 
by IStrubbe & Chiand (|2006|) with a more simplified approach. 



' However, these measurements include the so-called disc warp 
which, according to Golimowski et al. (2006), might be due to a blend 
of two separate, intrinsically thinner disc components inclined with re- 
spect to each other by a few degrees 

Our fit cannot give more quantitative information on the (e) and (i) 
dispersions, except that they are above the ~ 0.01 threshold value for 
the dynamically cold mode. 




best fit (<e> = 0.0035) 
_ nominal case {<e> = 0.1) 
observed profile 



70 80 90 100 
distance to the Star (AU) 



Fig. 12. The radial deprojected surface bright ness profile for 
the HR 4796 disc. Observed data are tak en by [Schn eider et 
(1999^ and are displayed in Fig.6 of IWahhai et al., (.20051 
Numerical profiles have been degraded to the resolution of the 
NICMOS/HST images, i.e., 0.12"=8AU. The best numerical fit 
is obtained for ArBR=16AU and <e) - 0.0035. 



For (5 Pic, this result should be reinvestigated in more detail in 
a more complete study taking into account the numerous addi- 
tional features observed in this system, but we believe our result 
to be robust regarding the global shape of this system's outer 
profile. Within the observational error bars, the nonimal result 
characteristics are also in relatively good agreement with outer 
disc profiles derived for several other debris discs, i.e. HD53143 
or HD32297 (Tab[T]l, especially when taking into account that 
strong anisotropic scattering could lead to somehow steeper pro- 
files in r""* (Tab[3]l. 

For all these systems, no additional sculpting mechanism is 
needed in order to explain the main ring's edge and the profile 
beyond it. Of course, something had to sculpt the parent-body 
ring in the first place, but such a sculpting could have occuiTed 
in the past, when for instance gas was still present, and no longer 
be active today. Another alternative is that the birth-ring's outer 
edge corresponds to the natural radial limit at which large plan- 
etesimals can form in the protoplanetary nebulae. 

There is however a group of systems (see Table [TJ which 
have slopes in the -4.5 to -5 range (HD139664, HD107146, 
HD181327), or even very sharp ring-like systems with brutal lu- 
minosity decreases in r^^ or r"^, represented by Fomalhaut and 
HR4796, which strongly depart from our nominal profile. Does 
this mean that these system's outer edges cannot be produced by 
"natural" internal collisional evolution and that "something else" 
is acting as a sculpting agent? The main question is here to see if 
these systems could fall into one of the categories giving alterna- 
tive profiles presented in Sec |3.2| especially the self-sustaining 
dynamically cold case. As a representative example, we consider 
the specific case of HR4796A.0 
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5.3. fitting a sinarp edge ring: HR4796A 

As for most debris discs, there exists no fully reliable observa- 
tional constraints on the vertical height and therefore dynamical 
excitation of the HR 4796A disc. From model fitting of scattered 
light and IR images, lAugereau et al I (jT999) argue that the scale 
height at 70AU (the location of the main ring) has to be less than 
7-8AU, or 2{i} ~ (e) < 0.1. Howev e r, fro m image fitting and 
theoretical constraints, iKenvon et alj (11999 ) find that the most 
likely value is H ~ 0.5AU, i.e. 2(0 ~ <e) < 0.007. It is thus 
plausible that this system falls into the dynamically cold cate- 
gory. 

We numerically explore a series of dynamically cold 
HR4796-like systems, fixing the ring's center at 70AU and tak- 
ing the radial extent of the birth ring ArBR as well as (e) as free 
parameters. The obse rved profile comes from HST/NICMOS 
(Schneider et al ][T999l) and is displayed in Fig. 6 of IWahhai et al.l 
(2005). Fig[T2] shows the best fit obtained in our parameter 
exploration. It corresponds to a dynamically cold case with 
(e) - 0.0035 and a radial width of the birth ring ArBR=16AU, 
which is close to the ^ 13 - 17AU d erived from observations 
(ISchneider et al.ll999t[Schneideil2001 )R As a comparison, we 
also display a "nominal" profile obtained for (e) = 0.1. The lat- 
ter is fully incompatible with the observed profile over the 80- 
120AU domain by almost 2 magnitudes in brightness. The best 
fit (dynamically cold case) compares well with the data in the 70- 
1 lOAU region and departs from the observed S B(r) profile only 
between 110 and 120AU, where the numerical profile becomes 
too flat. It is difficult to say how significant this discrepancy in 
these outermost lOA U is, as there are no error bars given in the 
IWahhai et al.l (l2005h plot. It is possible that confusion from the 
sky background enters at these distances. 

We must thus remain careful but it seems that there is at least 
a possibility for our dynamically cold model to provide an expla- 
nation for sharp outer edge discs like HR4796A. The question of 
how realistic this explanation can be is another issue. We have 
seen in Sec l5.1.2l that the (/) < 0.01 condition, although probably 
not generic with respect to planet-formation scenarios, cannot, in 
most cases, be explicitly ruled out by observations. 

There may however be a prediction made by the cold-disc 
scenario which could be observationally checked, i.e., the un- 
derabundance of grains in the j um to sub-millimetre ra nge F1. 

For HR4796, for instance, ^ Augereau et al.l d 1999 1) have per- 
formed detailed fits of the SED as well as of thermal and scat- 
tered light images. Their best fit implied that most of the geomet- 
rical cross section was contained in grains close to the minimum 
value for their size distribution, i.e. 5„„„ ~ lOjum. However, this 
fit was obtained assuming an imposed size distribution in s^^-^, 
so that these results cannot be used to rule out the possibility of 
higher i„„„ values for alternate size distributions. IWahhai et alj 
(l2005 h performed similar fits with partly more recent data and 
found that the effective size for grains within the ring is ~ 50 fim. 
This would be in relatively good agreement with our dynami- 
cally cold case. However, this value was obtained assuming a 



' The imaging data o n Fomalhaut do no t cover regions beyond ~ 20 
AU from the outer edge dKalas et alj2005l) . which is not enough to pro- 
ceed to a reliable numerical fit 

'° Note that since the disc's width is resolved in the near-IR 
( ISchneider et ai]|1999l : ISchneidej|200ll) . it could in principle be taken 
as a fixed input in our model. However, we kept it as a free parameter 
to check the validity of our fit. 

' ' While this paper was being reviewed, iLiseau et al ] ( l2008h reported 
the detection of a dust ring around q Eri, which could possibly be very 
confined and depleted from < 100/jm grains 



single size for the dust population, an obvious limitation the au- 
thors are fully aware of as they acknowledge the need for "phys- 
ically meaningful size distributions". It is a general character- 
istics of such global fits to have too many free parameters to 
constrain the dust composition as well as size and spatial dis- 
tributions in an unambiguous nondegenerate way. They usually 
have to rely on some starting assumption regarding the size dis- 
tribution. In this respect, it would be interesting to test these best 
fit models with alternate size distributions, as obtained from nu- 
merical collision-evolution studies, but this is an issue which ex- 
ceeds by far the scope of the present work. 

There exists however yet unpublished scattered light data 
(Schneider et al., http://nicmosis. as. arizona.edu:8000) indicating 
that the ring is uniformly red from V to H. This should mean 
that the emission is dominated by grains in the sub-yL/m to jim 
range. Should this feature be confirmed, then the dynamically 
cold state could probably be ruled out. Note however that this 
puzzling feature would be very difficult to interpret for any dy- 
namical model of the HR4796A ring. Indeed, it is not easy to 
explain how the luminosity, and thus the geometrical cross sec- 
tion could be dominated by unbound grains which should leave 
the ring on dynamical timescalesPl 

6. Summary and Conclusion 

We numerically investigate if the diverse outer surface bright- 
ness profiles observed in debris discs can be explained by the 
"natural" collisional evolution of belts made of solid particles 
ranging from planetesimal to micron-sized dust grains. We con- 
sider an initial ring of parent bodies with a razor-sharp outer edge 
and quantitatively examine to what extent the steady collisional 
production of small, radiation-pressure-affected grains modifies 
the initially perfect ring-like structure. 

Co ncurring with the pioneering results of lStrubbe & Chiang 
(I2006h . our numerical explorations have shown that for most 
"reasonable" parameters (mass, ring width, dynamical excita- 
tion) of a collisionally evolving debris ring, the surface bright- 
ness outside the ring naturally tends towards a standard profile, 
with no sharp drop at the outer edge and a mid-plane surface 
brightness profile oc r"^^. We confirm this result by simple but 
robust analytical considerations. This nominal result is in good 
agreement with the luminosity profiles observed in the outer re- 
gions of several debris discs, including j3 Pictoris and AUMic. 

Some observed systems however exhibit steeper radial lu- 
minosity profiles. Our numerical exploration shows that sharper 
outer edges and profiles steeper than that of the nominal case can 
only be obtained for two "extreme" cases: 

- 1) For a system with very high dust content (typically 
Mdusi ^ lOMe), the radial optical depth is raised to near 
unity for small grains. Most of these high-/? grains cannot 
travel out of the birth ring without suffering a collision. This 
allows the outer edge to sharpen. However, systems with 
such high radial optical depth are expected to wear down 
with time because of strong collisional erosion. Moreover, 
the existence of such very dusty systems is not backed by 
observations. This case does not appear as a realistic option. 

- 2) A longer survival of the sharp outer edge is achieved for 
systems with normal dust content that are dynamically cold. 

This seems to be a recurrent issue which is not limited to HR4796A: 
vast populations of unbound particles could also dominate the scattered 
luminosity of HD141569A ( lAugerea u & Papaloizou 2004), and might 
be present in the outer /3 Pictoris disc (lAugereau et alj200lT) 
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with typically (e) - 2{i) < 0.01. In this case, small grains 
are destroyed much more efficiently than they are created, 
leading to a depletion of this population. The system's opti- 
cal depth and luminosity are then dominated by large grains 
which do not leave the main birth ring, leading to a sharp 
outer edge. Even if this case might not correspond to the 
most generic debris disc configuration, it cannot be ruled out 
by observations and is thus a possible explanation to some of 
the observed systems. 

To numerically investigate the applicability of the dynami- 
cally cold case to real sharp-edge systems, we consider the spe- 
cific case of HR4796A. We find a reasonably good fit of this 
system's outer region luminosity profile with a dynamical ex- 
citation (e) ~ 0.0035. There is thus the possibility that such a 
sharp outer edge could be explained by the natural collisional 
evolution of a confined disc of large parent bodies. 

Acknowledgements. PT wishes to thank Jean-Charles Augereau and Alexander 
Kiivov for fruitful discussions. 
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